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ABSTRACT 

The prediction of torsional buckling of stiffeners is a necessary 
design skill in the planning of any framed structure or vehicle. The 
likelihood of torsional failure is increased when the stiffener is 
simultaneously loaded from perpendicular directions, a common loading 
condition for the frames of a ship at sea. This thesis examines first 
the case of a single lateral load on I beam and flat bar stiffeners 
and then the instance of combined lateral and axial loadings. An 
energy method is used and formulated are expressions for the critical 
buckling load for a single laterally applied force and a formula 
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NOTATIONS 



a Length of stiffener 

b Breadth of I beam flange 

B Angle of rotation of stiffener in the YZ plane 

Bo Amplitude of cosine function for B 

c Height of stiffener shear center above toe 

Cw Warping Constant 

d Additional lowering of the load P due to the load 

being applied a distance h/2 above center 
of the web 

E Young's Modulus = 2.068 x 10 8 KN/m 

G Shear Modulus = E/2.6 

h Height of stiffener 

Iz Moment of inertia of stiffener about web plane 

j,k,i Local rectangular coordinates at any point x 

on a stiffener 

J St. Venant's torsion constant 

L Dimensionless parameter 

M£ Moment about the l axis 

P Lateral load 

Per Critical lateral buckling load 

Q Axial end load 

Qe Euler buckling load 
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tf 



Flange thickness 
Web thickness 



tw 

u,v,w Displacements in the x,y, and z directions 
U Strain energy 

W Work 

x,y,z Rectangular coordinates originating at 
centroid of stiffener 
If Numerical factor 

C Membrane strain 

v Poisson's Ratio 
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INTRODUCTION 



As is the case with most large engineering structures, a ship's 
hull and its interior dividing bulkheads and decks are basically 
framework patterns of stiffeners covered by an outer shell. It is 
principally this lattice of stiffeners which determines the strength 
and structural stability of a vessel. Structural soundness and 
integrity are particularly critical for a ship because this 
engineering structure must operate in what can be the harshest natural 
environment on earth. Very seldom are the frames of a ship's hull 
subject to loading in a single plane; rather, it is most often a 
combination of loads caused by the interaction of the various natural 
and man-made forces from different directions which determine the 
total stress on a member. The subject of this thesis is to develop an 
expression relating the interaction of simultaneous lateral and axial 
Loading of a stiffener. 

Two of the more commonly found shipboard stiffeners, the I beam 
ind flat bar type, will be examined first under the action of a 
centrally applied lateral load and then under the combined actions of 
:he central lateral load and a per pind icular ly applied end axial load. 
Concentrated point loads were chosen because they represent the most 
.imiting and potentially dangerous cases. An energy approach will be 
ised in all instances. First the strain energy of the system will be 
ound and then the work done by the loading forces. By conservation 
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of energy, the variation of the strain energy will be equated to the 
work done. From the resulting expression will be developed, in the 
case of the single lateral load, a formula for the critical buckling 
load and, in the case of combined loading, a formula relating the 
interaction of the lateral and axial forces. 



It is hoped 
under more realis 
formulation and u 
military ship des 
lighter would cer 
the ultimate cont 
method of predict 
load may be susta 



that study into the problem of stiffener failure 
tic loading conditions such as this may aid in the 
pgrading of design codes used in civilian and 
ign, and any efforts that could make a ship safer or 
tainly be worthwhile. To this end it is hoped that 
ribution of this thesis is a reasonably accurate 
ing, given a value of axial load, how much lateral 
ined by a stiffener without buckling, or vice versa. 
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CHAPTER 1 



DESCRIPTION OF MODELS 

The two types of strengthening members that will be considered 
are the I beam type and flat bar type, to be known respectively as 
Model I and Model II. The coordinate system to be employed with both 
models is shown in figure (1-1). 



Z 




Figure (1-1) Coordinate System 



The origin of the 
of the beam, with u,v, 



coordinate system will be taken at the center 
and w being the displacements along the x, y, 
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and z axes. The deflections will be positive in the direction 
axes shown in figure (1-1). 

Torsional buckling is that failure mode characterized by a 
deflection and twisting of the stiffener in the YZ plane. This 
of deformation is shown in figure (1-2). 

Y 



X 

-V 
4 _ 




z 




Figure (1-2) Characteristics of Torsional Buckling 



of the 



type 
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1.1 Model I 



Model I as stated has been chosen to be a simply supported I 
beam. The case of a beam with symmetric flanges will be considered, 
and a lateral cross section of Model I is shown in figure (1-3). 



~^T 

_sk. 



T~ 





*1 




|< b 5 




h 



Figure (1-3) Model I Lateral Cross Section 



The properties of the Model I cross section are listed below: 
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c - Height of shear center above toe (Flange 

thickness t,f considered negligible compared 
to web height h) 

h 

c = 2 



Cw - Warping constant 



Cw = 




b 



3 



24 



J - St. Venant's torsion constant 



J = — -- ( 2bt f 3 + htw 3 ) 



Iz - Moment of inertia about the web plane 
(contribution of top and bottom 
flanges only) 



Iz 



b 3 t , 



2 ( 



12 



) 
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1.2. Model II 



Model II is the flat bar type stiffener as shown in figures (1-4) 
and (1-5). 




Figure (1-4) Model II - Flat Bar Stiffener 





Figure (1-5) Model II Lateral Cross Section 
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The properties of the flat bar cross section are as follows 



c - Height of shear center above toe 



c = 




Cw - Warping Constant 



Cw = 0 



J - St. Venant's torsion constant 



htw~^ 

3 



I z - Moment of inertia about the web plane 



Iz 



htw^ 

12 
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CHAPTER 2 



Model I - TORSIONAL BUCKLING UNDER POINT LATERAL LOADING 



2.1 Development of the Strain Energy Equation 

The case of a simply supported I beam of length a subjected to a 
point lateral load P in the center is illustrated in figure (2-1) 



P 




Figure (2-1) I Beam Under Point Lateral Load 



The general equation for the strain energy of a beam of length a under 
bending and torsion is given by (refernce 1): 



U = 




( -i^) 2 + GJ t- dB 



dx 



dx 



. 2 „ „ , d 2 B . 2 

■) + E Cw ( * ) 



dx 



□ 



dx 
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Considering only that portion of the beam to the right of an arbitrary 
point x, figure (2-2), the forces acting will be a single one P/2 in 
magnitude acting at x = a/2. If a local coordinate system with axes 
j, fc, H is assigned to the point x, the force P/2 will produce a 
moment about the & axis equal to: 



P 

M£ = 2 — ( a/2 - x ) 3 (2-1) 




with force P/2 acting to the right a distance a/2 - x 

Where B again is the angle of twist indicated in figure (1-2). It 
should be noted that B is assumed to be sufficiently small so that 
small angle approximations for sine and cosine are applicable. 



d 2 v 

The expression for then is (reference 1) 

dx 



d 2 v 

dx 2 



M 

& 

El 
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For B the simple choice B = Bo cos ~ will be made. The cosine 
expression is selected because it is assumed that the stiffener will 
be mounted such that it is restrained from twisting at the ends x = 
a/2 and x = -a/2. This is a reasonable assumption in the case of 

stiffeners fixed between decks or frames of a ship. 



d 2 v 

Utilyzing the appropriate terms for M 0 and B, 7 — becomes: 

* dx 2 



d 2 v 



P , a . _ 7TX 

5 - ( ■= - x Bo cos 

j z 2EI 2 a 

dx z 



( 2 - 2 ) 



Inserting this expression into the formula for strain energy U yields: 



U = 



- a/2 p 2 Bo 2 



2 ttx , a ,2 , __,Bott.2 • 2 ttx 

JET - COS — < 2 - X) + GJ( — > Sln T 



, _ 2 , it . 4 2 ttx 

+ ECw Bo ( - ) cos 

a a 



dx 



Carrying out the integration it 

2 2 2 

U = P ■ B ° — a 3 ( - + ) + GJ ( 

El 192^ 

z 

2.2 Development of the Virtual 



becomes 



Bott ,2a , „„ „ 2 

) t + ECw Bo 

a 4 



Work Equation 



tt. 4 a 
a 4 



Work is defined as the movement of a force through a distance. 
In the case of Model I then the virtual work done is obtained by 
multiplying the force P times the summation of all the components of 
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motion in the 2 direction for elements of the beam between x = 0 and x 



= a/2. It is not necessary to sum vertical components over the entire 
length of the beam, x = -a/2 to x = a/2 , because each half moves the 
same total distance. The expression for the vertical distance moved 
by the centroid of a cross section of the beam at the arbitrary point 
X is (reference 1): 



This equation would suffice alone if the force P were applied at the 
centroid of the cross section at X = 0 ; but since we are considering 
the more realistic case where the load is applied at the outside 
surface of the beam, a distance h/2 above the centroid, the expression 
must be slightly ammended. Due to the twist of the section at X=0, 
the load will drop through an additional distance d = h/2(l - cos B*, 0 ) 
as shown in figure (2-3). 





and the work equation then becomes: 



a/2 



W = P 




^ - x ) dx 
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p 




Figure (2-3) Additional lowering of P due to load 
being applied a distance h/2 above centroid 



By the half angle relation, 



2 

1 - cos B xco = 2 sin B^^/2 



= 2(W2) 2 

= B?.o/2 



Therefore d becomes approximately: 

d - (h/2 ) (bL/2) 

= (h/4 ) Bo 2 cos 2 0 
= (h/4 ) Bo 2 
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This term is now multiplied by P and added to the work equation to 

give a more accurate expression for the work done. 
a/2 



W = P 



5 . 



. 2 

n d v , a > , 

B — ( j — x) dx + 

dx z 



Ph D 2 
-t-Bo 



(2-3) 



a/2 



= P 



Bo cos 



7TX 



2 

, d v > , a >, ,Ph 0 2 

( tt ) (-s — x) dx + Bo 

dx 4 



2 

d v _ P 

Again — 2 2EI~ 

dx z 



(~ 2 — - x ) B , and 



when applied the work equation becomes: 



W = 



2 2 
P Bo 

2EI 



a/ 2 



cos 



TTX 

a 



N 2 , ^ Ph _ 2 

- x ) dx + — Bo 



Integrating this expression gives: 



r7 P 2 Bo 2 3 , 6 + it 2 , , Ph ^_2 

W = a ( — ) + — r— Bo 



El 



96 



2.3 Determination of the Critical Buckling Load 



To determine the value of P that will cause failure of the 
stiffener, the critical buckling load Per, the principle of 
conservation of energy will be employed. This concept states that 
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when a member is gradually loaded, the kinetic energy is zero and the 
work done by external forces is equal to the strain energy (reference 
2 ) : 

6W = 6U 



To determine 5W and 6U calculus of variations with respect to Bo is 
applied to the work and strain energy equations. First to the work 
equation : 



W ( Bo + 6Bo) 



= (Bo + 6 Bo) 



2 3 2 

Pa . 6 + ir 

EI qc 2 

z 96 tt 



= W ( Bo ) + 2 Bo 6 Bo 



B 



P 2 a 3 , 6 + * 2 . 
EI z 96 tt 2 




5W 



2 Bo 6 Bo 



r e!Z 

EI 




And in like manner to the strain energy equation: 
U (Bo + 5Bo) = (Bo+-5Bo) 



2 3 2 ~ 

2 | P a / 6 + ^ \ , r* -r \ ^ ^ i r-» yv , ( \ ^ ^ 

IT ( r7nT )+GJ( a' 4 + E c w( — ) 5 



= U (Bo) +B06B0 



D 

0 



z 192 tt 



<^7T 2)+ 75 2(gj + ECW( 



6U=2Bo6Bo 



z 192 tt' 
.2 



a'L] 



|y^-( 6+TT - y - )+ JZ (GJ + ECw ( 
EI Z 1927T 2 4a 






The two variational expressions may now be equated, and the formula 
for the critical buckling load derived. 

6W = 6U 

24 



2Bo6Bo 



2 3 2 

Per a . 6+tt Perh 

' -1 ’ + A 



El 



96tt‘ 



— )+ ra (GJ + EC «<i> 2 ] 

z 192 tt — 1 



z 

Per^a^ , 6 +tt^ n , Per h tt ^ ^ ,1,2 

( -3) + — 4^ |^GJ -f ECw ( — ) 



El 



192tt‘ 



3 - 



This quadratic may be solved for Per, and the result expressed in the 



form : 



Per = 1 - ( 



24 



6+tt‘ 



W- 



EI . 

— 2 + l 

GJ 2 



. 48 tt^ 2 h. 2 EI z 

2 a GJ 

6 + tt 



+ ) ( 1 + ECW 



6+tt 



GJ 



( l ) 2 ) 

a. 



1/2 J y'EIGJ 



(2-4) 



This formula corresponds directly with equation (6-18) of reference 1 
for similar geometry and loading which, in the notation of this 
thesis, is: 

y^i gj 

Per = ~ 



The expression within braces in equation (2-4) equates to the non- 
dimensional factor of (2-5). Figure (2-4) shows the close agreement 
between the values of t calculated from equation (2-4) with those 
listed in tabular form in Table 6-5 of reference 1. 
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Figure (2-4) Agreement of calculated values of Y from eq (2-4) 
with tabular values listed in reference (1) 
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CHAPTER 3 



MODEL II - TORSIONAL BUCKLIN3 UNDER POINT LATERAL LOADING 



3.1 Development of the Strain Enerqy Equation 



Model II again was the flat bar type stiffener illustrated in 
figure (1-4). In this chapter the same type of loading situation as 
that of chapter 2, a point load coincident with the Z axis and applied 
at the upper, surface of the stiffener, will be examined. A pictorial 
representation of the loading will appear identical to that shown in 
figure (2 -1 ) . 

The equation for the strain energy of a flat bar in lateral 
bending and torsion can be obtained in exactly the same manner as that 
for the I beam. The only difference is that the warping riqiditv 
term, Cw, vanishes. This is because the equation for Cw of an I beam 



and the flat bar stiffener can be likened to an I beam with the width 



b of the top and bottom flanges equal to zero. The strain energy 
expression then is: 

a/2 t — ? ? 2 



U = 



(S ^2> 2 + GJ '-if-' 






dx 



The simple choice for B, Bo cos - — , is again made and the 

d 2 3 

expression for 2— X is then identical to (2-1): 

dx 2 



d 2 v 



dx 2EI 



P , a . „ TTX 

— ( - x) Bo cos — 

A 3 



The complete strain energy eauation is: 



a/2 



U = 



r* 2 2 2 

P Bo COS 

4EI a 

— z 



( - 
V 2 



. 2 , _ T , Bott> 2 . 2 

x) + GJ ( ) sm 



d 



dx 



which when integrated becomes: 



2 2 2 2 
P Bo 3 , 6 + tt , ,Bott, a 

U = a ( ) + GJ ( ) — 



4EI 



192tt 
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3.2 Development of the Virtual Work Equation 



The work equation is found again by defining the work done as the 
product of the load P times the summation of all vertical components 
of motion over one half of the bar. The expression for the vertical 
distance moved by the centroid at any point along the X axis is 
identical to the I beam case: 



B 



d 2 v 



d x' 




x) 



Also, a correction term must once more be added to the work equation 

to account for the load being applied at the upper surface of the 

plate rather than at the centroid. The ultimate result is that the 

work equations for model I and model II are identical: 

a/2 

T7 P 2 Bo 2 C 2 ttx ,a . 2 , , Ph _ 2 

W = \ cos (— - x) dx + — Bo 

2EI \ a 2 4 

Z J 

v> o 

P 2 Bo 2 _ 3 ,6 + it 2 , , Ph „^2 

= a ( 2 — ) + — Bo 

El 96tr 4 

z 

3.3 Determination of the Critical Buckling Load 



By setting the variation of the work done equal to the variation 
of the strain energy the value of Per is once again determined. Since 
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the work expression for the flat plate is identical to that of the I 
beam, so too is the equation for the variation of the work: 



2 3 2 

rn On rr, | Pa ,6 + I , , Ph 

6W = 2Bo<5Bo | ( ) + — 

El 96 tt 4 

z 



[ 



1 



While that for 6U is slightly different owing to the abscence of Cw 



U ( Bo + 6 Bo) = (Bo + 6 Bo) 



(§.. + + G J (V ^ 

EI 192tt a 4 



J 



[ 



6U = 2 Bo 6 Bo | P a ( 6 + - j ) + GJ (-) | 

El 192tt a 4 

z 



] 



Equating the work and strain energy and solving for Per yields: 



2Bo6Bo 



6 W = 6U 

2 3 2 ~ 

Per a ,6 + 7T ^ , Ph 

k 2 ’ + 

El 96 tt 4 



= 2BooBo 



2 3 2 2 

Per a ,6 + ir . . ,tt. a 

~~ ( " o ~ 2 ) + GJ a 4 

El 192 tt 



( 6 + tt 2 } _ £h + if. GJ = 0 

El 192tt 4 4a 

z 



When solved this quadratic gives for Per an expression which may be 
written in a form similar to equation (2-2): 



30 



Per = 



{ 



* ys - + 1 

6 + 7T a V GJ 2 



“ , Q 2 2 ,2 El, 

(il) — 



6 + tt‘ 



GJ 



+ 



4 

( 192tt 

6 + TT 




1/2 




1 / El GJ 

* z 

2 

a 



(3-1) 



TOO TT" 

For typical stiffeners the quantity ( ) will be several 

2 6 + ei 

orders of magnitude greater than ( — 4_8_rr — ^ ^h j2 — _ — z . Therefore, 

6 + TT Z a GJ 

purposes of a design formula it is sufficiently accurate to ignore 
this latter quantity and rewrite equation (3-1) as: 



Per 



1 / 192tt 4 ,1/2 

' 2 ' 

2 6 + TT 



24- 



6 + 



r > * 





EI GJ 
z 

2 

a 



(3-2) 



The formula for Per for a flat plate stiffener subjected to a 
point load in the center of the span is given in reference (1) by 
equation (6-37), which in this thesis' notation is: 



Per = 16.94 






EI z GJ 



(1 





) 



Now to two-decimal accuracy equation (3-2) may be rewritten as 



Per = 17.16 



sf. 



EI Z GJ h 

^2 — (1 " - 87 a 



l 



EI 



GJ 



The formulas are in agreement within 1.3%. 



for 
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CHAPTER 4 



MODELS I AND II - TORSIONAL BUCKLItG UNDER COMBINED 
LATERAL AND AXIAL POINT LOAD INS 



4.1 Description of Loading Condition 



In chapters 2 and 3 the failure by torsional buckling of a beam 
subject to a lateral load only has been studied; now suppose that an 
additional axial load is added at the ends of the beam as shown in 
f igure (4 -1 ) . 



P 




Figure (4-1) Combined Lateral and Axial Loading of a 
Stiffener 
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The stiffener is again simply supported and laterally loaded by a 
force P applied at the upper surface in the center of the span, but an 
additional force Q acts at the centroid of each end to compress the 
stiffener inward. The two forces are acting at right angles to each 
>ther; P along the Z axis and Q along the X axis. This dual axis type 
oading would seem to be much more representative of the kind of 
oading conditions to be encountered by the frames of a ship. For 
xample, two commonly encountered conditions by a ship at sea in large 
aves are hog and sag, shown in figure (4-2). 





Figure (4-2) Hogging and Sagging Congitions in 
Large Waves 
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Here the ship is supported by the crest of a wave either amidships 
(hogging) or at the ends (sagging). In both cases the vertical 
transverse frames of the ship, a portion of which are shown as dashed 
Lines in figure (4-2), will be placed in compression, i.e. a Q force 
applied. At the same time then, the sides of the ship will very 
Likely be receiving a lateral force P from other waves in a confused 
sea or in a worst case from collision with another ship, floating 
abject, or possibly from within the ship from unsecured cargo 
auf feting about. How do the two forces interact, and what effect does 
the addition of an axial force have upon the capability of the 
stiffener to accept a lateral load? These are the questions this 
:hesis will attempt to answer utilizing the strain energy approach 
/hich was used in the previous chapters for a lateral load alone. 



1.2 Development of the Strain Energy Equation 



In solving for the strain energy of a stiffener subjected to both 
.ateral and axial loads the case of model I, an I beam, will be 
:onsidered first; so the energy equation will contain a warping 




EI z '^%> 2 + GJ ( -£ r > 2 + ec w <^! )2 



dx 



dx 



dx 
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Again for B, the angle of twist, the choice Bo cos _JS is made. The 

2 a 

d v 

expression for will however, no longer be the same. Taking again 

dx^ 

an arbitrary point x as shown in figure (2-2), there will still be 
present to the right of x the force P/2 creating a moment about the 
local £ axis equal to the expression in equation (2-1): 



M£ = | " x ) B 



( 2 - 1 ) 



There is now though, another force to the riqht of x adding to this 
moment about?. , the axial force Q. Looking down upon the stiffener 
from above it is seen that this addition to M i is equal to the force Q 
times the lever arm -v as shown in figure (4-3): 




Figure (4-3) Additional moment Qv provided by 
axial loading 



The total moment M £ then is now: 
M£ = | (| - x) B + Qv 
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/ 



ind the expression for 



becomes : 



d 

dx 




Mil 

El 

z 



P 

2EI 

z 



a 

2 



X ) Bo cos 



TT X 

a 




(4-1) 



tfe now have a linear second order differential equation which must be 
solved for v in order that the strain energy may be found. We are 
Looking for a particular solution to this differential equation, and 
:hus it will be of the form (reference 3): 



Vp = C-^ cos 



TT X 



+ sin 



TT X 



+ C 3 x 



COS 



TT X 



C. x 
4 



sin 



TTX 



(4-2) 



) if f erent iat ing this expression for v twice with respect to x yields: 

2 h 

d Vp TT TT . TTX TT , TT _ . . TTX 

— x = - (2C. - - C, cos - - (2C 0 + - C 0 ) sin 

,2 a 4 a 1 a a 3a2 a 

dx 



.TT 2 „ TTX , TT . 2 _ .TTX 

<i> C 3 X C ° S — ' ( 5> C 4 X Sin T" 



(4-3: 



[f v is a solution to equation (4-1) the following is true: 
P 



d Vp P ,a * „ TTX , 0 / ^ TTX , „ . TTX 

— (■= - x) Bo cos 1 -=v (C. cos + C_ sm 

, 2 2EI 2 a El 1 a 2 a 

dx z z 



„ TTX „ . TTX 

+ C x cos + C. x sin ) 

3 a 4 a 



(4-4) 
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By equating coefficients of the trigonometric functions in equations 
(4-3) and (4-4), four equations can be found to be solved for the four 



constants through C^. They are: 



C 1 „ 2 



-a PBo 



4a Q + 4EI t 
z 



C 2 = 



-a El tt PB o 
z 



2 2 2 

(a^Q + EI z *-V 



c = 

3 .2 



a 2 PBo 



2a Q + 2EI tt ' 
z 



c 4 = 0 



Inserting these constants into equations (4-2) and (4-4) gives for v^ 
and its second derivative: 



V P = ( — 



-a PBo 



4a Q + 4EI tt 
z 



, tt x , a El ttPBo . . tt x 

tt ) cos — - ( z ) sm — 

2 a 5 2 2 a 

(a Q + El tt ) 
z 



+ ( 



a 2 PBo 



2a Q + 2EI tt 
z 



. TT X 

■S' ) X COS 

2 a 



(4-5i 



and 



,2 

d v 



dx‘ 



P , a . _ ttx Q 

= 2EI < 2 - X) 20 005 T + ET 



-a 3 PBo 



4a Q + 4EI rr 
z 



, TT X 

TT )cos — 
2 a 



a El ttPBo 

/ z . . ttx , , 

— 5 TT sm — — + (- 



2 2 2 
(a Q + El tt ) 
z 



2 

a PBo . tt x 

* )x cos 

z a 



2a Q + 2 EI z tt 



(4-6) 



,2 

d v 



dx 



PBO TTX 

"ei cos ~ 



a 3 PBoQ 



( - - 

1 4 .2 



a 3 Q 



,) + (- 



4a"Q + 4EItt 2 2a 2 Q + 2EI tt 

z 



B 




(a Q + El tt ) 

v Z 



TT X 

277 sin — ~ 



(4-7) 



For simplicity the coefficients in equation (4-7) shall be labeled 
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, B , and C as indicated, and (4-7) then becomes: 



d 2 v 

dx 2 



PBo 



El 



T fY 7T V 

cos — — (A + Bx) - C sin 



nd 



( d - v ) 2 = 
dx 2 



, PBo ,2 2 ttx , ,2 , , „2 2 N 

( — ) cos ( A + 2ABx + B x ) 

LX 3 



- 2 — XjjL — (AC + BCx) sin — — cos 2 — — — + C 2 sin 2 — — 
El a a a 



(4-8) 



ow with equation (4-8) there is sufficient information to solve the 



train energy equation. 



U =' 



-a /2 

J o 
-a /2 



El ( — 5 ) 2 + GJ ( ) 2 + ECw ( -X- B ' 2 

2 ^ Cl X ci “ 



') 



dx 



dx 



dx 



r 2 n 2 

p Bo 



El 



2 "x / ,2 . 2 2. 

cos ( A + 2AB x + B x ) 



- 2PBo ( AC + BCx ) sin cos 



TTX 



„2 . 2 ~x , , Bo r 1 2 . c x 

+ El C sin + GJ ( •) sin 

z a a a 



„ _ _ 2 , 7T . 4 2 TTX 

+ ECwBo ( - ) cos 

a a 



dx 



(4-9) 



arrying out the integration of equation (4-9) and reinserting the 
roper quantities for A,B, and C gives as a final expression for 
Crain energy: 
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